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ABSTRACT 

We conducted a comprehensive study on the shell structure of the Cygnus 
Loop using 41 observation data obtained by the Suzaku and the XMM-Newton 
satellites. To investigate the detailed plasma structure of the Cygnus Loop, we 
divided our fields of view into 1042 box regions. Prom the spectral analysis, the 
spectra obtained from the limb of the Loop are well fitted by the single-component 
non-equilibrium ionization plasma model. On the other hand, the spectra ob- 
tained from the inner regions are well fitted by the two-component model. As a 
result, we confirmed that the low-temperature and the high-temperature compo- 
nents originated from the surrounding interstellar matter (ISM) and the ejecta 
of the Loop, respectively. 

Prom the best-fit results, we showed a flux distribution of the ISM component. 
The distribution clearly shows the limb-brightening structure, and we found out 
some low-flux regions. Among them, the south blowout region has the lowest flux. 
We also found other large low-flux regions at slightly west and the northeast from 
the center. We estimated the former thin shell region to be ~ 1°.3 in diameter 
and concluded that there exists a blowout along the line of sight in addition to the 
south blowout. We also calculated the emission measure distribution of the ISM 
component and showed that the Cygnus Loop is far from the result obtained by 
a simple Sedov evolution model. From the results, we support that the Cygnus 
Loop originated from a cavity explosion. The emission measure distribution also 
suggests that the cavity-wall density is higher in the northeast than that in the 
southwest. These results suggest that the thickness of the cavity wall surrounding 
the Cygnus Loop is not uniform. 

Subject headings: ISM: abundances — ISM: individual (Cygnus Loop) — supernova 
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remnants — X-rays: ISM 
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Introduction 



The supernova (SN) explosions blow out the various heavy elements generated by the 
nucleosynthesis process inside the progenitor stars. Meanwhile, the blast wave generated by 
the SN explosion sweeps up and heats the interstellar matter (ISM). It forms a characteristic 
structure of each supernova remnant (SNR). In this way, the morphology of the shell 
structure provides information about the ambient density of the ISM. 

The Cygnus Loop is one of the brightest SNRs in the X-ray sky. Its age is estimated 
to be ~ 10, OOOyr and the distance is comparatively close to us (540pc; Blair et al. 
2005). The large apparent size (2.5°x3.5°; Levenson et al. 1997) enables us to study the 
plasma str ucture of the Loop. The origin of the 



explo sion (IMcCray fc Snow 



I979I : 



Hester fc Cox 



1986 



Hester et al. 


1994 




Levenson et al. 


star ( 


Levenson et al. 


1998 


) and 



some X-ray studies also estimated the progenitor mass to be 12-15M0 {e.g., Tsunemi et 
al. 2007). From the morphological point of view, the Cygnus Loop is a typical shell-like 
SNR and it is almost circul ar in shape. However, a la rge break is seen in the south. 



known as "blowout" region (lAschenbach fc 



not well understood. 



Aschenbach fc Leahyl ( 



Leahv 



19991 ). The origin of the "blowout" is 



I999I ) explained this extended structure as 



a breakout into a low er density ISM. On the other hand, based on a radio observation. 



Uyaniker et al. 



(I2OO2I ) suggested the existence of a secondary SNR in the south. Some 



other rad io observation s supp ort this conclusion (jUyaniker et al. 



Recently, 



Uchida et al. 



2004 



Sun et al. 



20061 ). 



(I2OO8I ) observed this region with the XMM-Newton observatory and 
found that the X-ray spectra of this region consist of two plasma components with different 
electron temperatures. Judging from the plasma structures and the metal distributions, 
they concluded that the X-ray emission is consistent with a Cygnus Loop origin and that 
each plasma component is derived from the ejecta and the cavity material, respectively. 
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They also showed that the X-ray shell is thin in their fields of view (FOV) and concluded 



that the ori gin of the blowout can 
proposed by lAschenbach fc Leahyl (Il999l ) 



3e ex plained as a breakout into a lower density ISM as 



It is natural to c onsider that such a breakout may also exist along the line of sight. 



Tsunemi et al. 



(120071 ) observed the Cygnus Loop along the diameter with XMM-Newton 
and argued about it. They divided their FOV into a north path and a south path and found 
the thin shell region to be 5' in the south path and 20' in the north path. They estimated 



this thin shel 



Kimura et al. 



region to have a diameter of 1°, centering on a = 20^9™!!*^, 6 = 31°05'20". 



(120091 ) expanded their observation northward with Suzaku and found that the 



fiux of the swept-up m atter in southwest is about a third of that in northeast. The width 



of this region is ~ 50'. 



Kimura et al. 



( 120091 ) presumed that there is a blowout region in the 



direction of our line of sight in the middle west of the Loop. 

In this paper, we used 41 observation data obtained by the Suzaku and the XMM- 
Newton observatories. We reanalyzed all the data to conduct a comprehensive study on the 
shell structure of the Cygnus Loop. 



2. Observations 

We summarized the 41 observations in table [H All the data were taken between 2002 
and 2008. The observing regions are shown in figure [T] top. The circles and rectangles 
represent the FOV of the XMM-Newton MOS and the Suzaku XIS, respectively. 

All of the Suzaku data were analyzed with version 6.5 of the HEAsoft tools. For 
reduction of the Suzaku data, we used version 9 of the Suzaku Software. The calibration 
database (CALDB) used was updated in July 2008. We used revision 2.2 of the cleaned 
event data and combined 3x3 and 5x5 event files. We applied the spaced row charge 
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injection (SCI) method flPrigozhin et al.ll2008h to PI, P2, P3, P4, P5, P6, P7, P9, PIO, 
Pll, P18, P19, P20, P21, P22, P23, P24 and P25 data. This method reduces the effect of 
radiation damage of the XIS and recovers the energy resolution, for example, from 205±6eV 
to 157±4eV at 5.9 keV. In order to exclude the background flare events, we obtained good 
time intervals (GTIs) including only times at which the count rates were within ±2a of the 
mean count rates. 

Since the Cygnus Loop is a large diffuse source and our FOV are filled with the 
SNR's emission, we could not obtain background spectra from our FOV. We also had no 
background data from the neighborhood of the Cygnus Loop. We therefore applied the 
Lockman Hole data for background subtraction. We reviewed the effect of the galactic 

The flux of the GR XE at / = 62°, \b\ < 0°.4 is 6 x lO^^^g^g 



ridge X-ray emission (GRXE 
cm~^s~^deg~^ (0.7-2.0 keV) (ISugizaki et a. 
b = -8°. 6) is located outside the FOV of 



20011). Although th e Cygnus Loop (/ = 74°, 



Sugizaki et al. 



(1200 ll ). this value gives us an 



upper limit of the GRXE at the Cygnus Loop. Meanwhile, the flux of the Cygnus Loop 
is estimated to be 8.2 x 10~^°erg cm^^s^^deg^^ (0.7-2.0 keV), assuming that the Cygnus 
Loop is a circle with a diameter of 3°.0. Therefore, we concluded that the effect of the 
GRXE on the Cygnus Loop is vanishingly small. The solar wind charge exc hange (SWCX) 



is als o considered to correlate with the soft X-ray background below 1 keV ( iFujimoto et al 



20071 ). However, in terms of the Cygnus Loop, we consider that the SWCX is negligible 
because of the prominent surface brightness of the Loop. Thus, the Lockman Hole data 
obtained in 2006, 2007 and 2008 were applied for background subtraction. We selected data 
whose observation dates were close to those of the Cygnus Loop observations. Since there 
were no photons above 3.0 keV after background subtraction, the energy ranges of 0.2-3.0 



keV and 0.4-3.0 keV were used for XISl (back-i 



(front-illuminated CCD; FI CCD), respectively (jKoyama et al. 



luminated CCD; B 



2003). 



CCD) and XIS0,2,3 
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AU the XMM-Newton data were processed with version 7.1.0 of the XMM Science 
Analysis System (SAS). The current cahbration files (CCFs) used were updated on 2008 
June. We used data obtained with the EPIC MOS and pn cameras. These data were 
taken by using the medium filters and the prime full-window mode. We selected X-ray 
events corresponding to patterns 0-12 and fiag = for MOS 1 and 2, patterns 0-4 and 
fiag = for pn, respectively. In order to exclude background fiare events, we determined 
the GTIs in the same way as those of the Suzaku data. After filtering the data, they were 
vignetting-corrected by using the SAS task e yigweight. For backgro und subtraction, we 



Read fc Ponman 



( 120031 ) for a similar reason 



employed blank-sky observations prepared by 
with the Suzaku case. After the background subtraction, the energy range of 0.3-3.0 keV 
was used for each instrument. 



3. Spectral Analysis 

To investigate the plasma structure of the Cygnus Loop, we divided the entire FOV 
into small box regions. In order to equalize the statistics, we initially divided all images 
of XISl or M0S2 into two parts; if each divided region had more than 10,000 photons, it 
was once again divided. In this way, we obtained 1042 box regions. Each region contained 
5,000-10,000 photons for XISl and M0S2. The side length of each box ranges from 2.2' 
to 14'. Therefore box sizes are not smaller than the angular resolution capability of the 
Suzaku XIS. We grouped 1042 spectra into bins with a minimum of 20 counts such that 
statistics are appropriate. In order to generate a response matrix file (RMF) and an 



ancillary response file (ARE), we employed xisrmfgen (jlshisaki et al.ll2007l ) and xissimarfgen 



for the Suzaku data, rmfgen and arfgen for the XMM-Newton data. 

Firstly, we fitted all the spectra by a single-component variable abundance non- 
equilibrium ionization (VNEI) model. We employed TBabs (Tiibingen-Boulder ISM 
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absorption model; Wil ms et al. 2000) and VNEI (NEI ver.2.0; Borkowski et al. 2001) in 



XSPEC version 12.5.0 flArnaud 



19961 ). In this model, the abundances for C, N, O, Ne, 
Mg, Si and Fe were free, while we set the relative abundances for S to the solar value equal 
to that of Si, Ni equal to Fe. The other elements were fixed to their solar values. Other 
parameters were all free such as the electron temperature, kTg, the ionization timescale, 
T (a product of the electron density and the elapsed time after the shock heating), and 
the emission measure, EM (= J rienndl, where Ug and tt-h are the number densities of 
electron and hydrogen and dl is the plasma depth). We also set the absorption column 
density, A^h, to be free. As a result, the spectra from the limb regions are well fitted by 



the single-component V^ 
and the southeast limb ( 


EI model. As shown by earlier observations of the northeast 


Tsunemi et al. 


2007; 


Kimura et al. 


2009; 


Uchida et al. 


2009a: 



Tsunemi et al. 



20091 ). the spectra obtained from the limb regions of the Cygnus Loop are 



typically described by a single-component VNEI model. 

On the other hand, the spectra from the inner regions are generally not fitted by 
the single-component VNEI model. From earlier obs e rvatio ns of the northeast to the 



southwest regions along the diameter. 



Tsunemi et al. 



(120071 ) found that the spectra from 



the inner regions of the Cygnus Loop consist of a two-component VNEI plasma. They 
concluded the plasma structure of the Cygnus Loop as follows: the high-ZcT p ejecta 



Uchida et al 



fl2009bh showed that 



component is surrounded by a low-fcTg ISM component, 
the two-component VNEI model is wholly applicable to the inner regions of the Cygnus 
Loop. Therefore, we next intended to give an additional high-/cTe VNEI component to 
the single-component VNEI model. In this model, we fixed the metal abundan c es of 



the low-A;Te component to the values obtained from the result of 



Tsunemi et al. 



(120071 ) 



since the model whose a 
results 



Tsunemi et al 



jund ances set all free could not obtain the physically meaningful 



( 120071 ) showed the relative abundances to the solar values of the 



ISM component as follows: C=0.27, N=0.10, 0=0.11, Ne=0.21, Mg=0.17, Si=0.34, S=0.17, 
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Ni)=0.20. In addition, they fixed other elements to the solar values (lAnders &: Grevesse 



19891 ). Meanwhile, in the high-A;Te component, the abundances for O, Ne, Mg, Si, and Fe 
were free, while we set the abundances for C and N equal to O, S equal to Si, Ni equal to Fe. 
Other elements were fixed to their solar values. The other parameters such as fcTe, r, EM, 
and A^H were all free. We applied both single-component VNEI model and two-component 
VNEI model to all the spectra and determined which model is acceptable by using the 
F-test with a significance level of 99%. As a result, roughly < 0.80i?s of the northeast region 
and < 0.85i?s of the southwest region need an additional component, where Rs is a shock 
radius. Here, we define the "limb observations" as the regions where the single- component 
VNEI model is acceptable and the "inside observations" as the remaining regions. 

Figure [2] shows two example XISl spectra. The spectral extracted regions are shown in 
figure [H Both regions are located at the inside observations and the two-component VNEI 
model is applicable. The bottom two panels show the best-fit results of the two-component 
VNEI model. Blue and red lines represent the high-/cTe and the low-fcTg component. We 
also show the result with the single-component VNEI model at the top two panels for 
comparison. The best-fit parameters are shown in table [2l These results show that the 
reduced values are significantly improved with the two-component VNEI models. 



4. Discussion 

4.1. Temperature distribution of the Iow-A;re component 

All the spectra are well fitted by either the single-component VNEI model or the 
two-component VNEI model. From the best-fit parameters of the inside observations, we 
found that the electron temperature of the low-kTe component is almost uniform. The 
averaged value is 0.23 keV {a = 0.08 keV) and it is sufficiently lower than that of the 
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high-/cTe component (0.52 keV, a = 0.17 keV). The temperature of the low-kTe component 
is close to that of the hmb observations (0.29 keV, a = 0.07 keV). Therefore we collectively 
call these components "low-ZcTg component" hereafter. 

Figure [3] shows our FOV and the electron temperature distribution of the low-fcTg 
component overlaid with the white contour from the ROSAT HRI image. The averaged 
value is ~ 0.28 keV and it ranges from 0.12 keV to 0.35 keV. Meanwhile, the temperature 



of the high-fcTe compo nent ranges from 0. 



previous observation s (Tsunemi et al 



2007 



^eV to 0.9 keV, which is consistent with t 



he 



Katsuda et al. 



Uchida et al. 



2009 



clearly separated. 



2008b : 



Kimura et al.l 12009 



a). Then, we confirmed that the temperature of each component is 



Uchida et al. 



(j2009bl ) also showed that the temperature distribution of 
the high-fcTg component is not uniform and that it is lower in the southwest part than 
that in the northeast part. On the other hand, the temperature of the low-A;re component 
is relatively uniform (see figure [3]). The detailed distribution shows the temperature near 
the center is lower than that of the surroundings. We also found that the temperature 
distribution is seamless at the boundary between the limb observations and the inside 
observations. Therefore, the low-A;Te components of these regions must have the same 
origin. The spectra from the limb observations are obviously swept-up ISM origin, and 
thus, we concluded that any low-ZcTg component originates from the ISM component. 



4.2. Line-of-sight Shell Structure of the Cygnus Loop 

Taking into account the age of the Cygnus Loop the reverse shocks should have already 
reached its center. Therefore, on the assumption that the density of the ejecta-origin plasma 
is homogeneous, the X-ray fiux depends exclusively on its plasma depth. In figure [21 the 
blue line represents the high-fcTe component of the two-component VNEI model. Since the 
region-A and the region-B are located at the same radial distance from the center {R ~ 50', 
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where we define i? as a distance from the "geometric center" determined by Levenson et 
al. 1998), they should have almost the same plasma depths. Accordingly, the fiuxes of 
the high-fcTg components are actually not so different, while the spectral extracted regions 
are separated. Meanwhile, the contributions of the low-A;Te components are quite different 
as shown with the red lines in figure O From the bottom left panel of figure O the fiux 
of the low-A;Te component in the region-A overwhelms that of the high-ZcTg component at 
0.2-1.0 keV. On the other hand, the contribution of the low-A;Te component in the region-B 
is clearly smaller than that in the region-A. Such a difference should be attributed to the 
difference of the surrounding shell of each region. The value of the flux is proportional 
to EM (oc n^/), which means that the surface brightness is sensitive to the change of the 
density and the plasma depth there. 

In order to estimate the ambient density of the Cygnus Loop, we calculated the fluxes of 
the low-kTf. component from all regions. FigurelUleft shows the 0.2-3.0 keV flux distribution 
of the low-A;Te component. We also show that of the high-fcTe component at the right panel. 
The flux distribution of the high-fcTg component is relatively uniform compared with that 
of the low-/cTe component. It reflects that the ejecta component uniformly flUed inside the 
Loop. In contrast, from the left panel, we clearly see the "limb-brightening" which reflects 
the spherical shell structure. Therefore, we conflrmed that the low-kTg component comes 
from the surrounding ISM. We also found that the northeast flux is higher than that in 
the southwest. It suggests that the density is higher in the direction of the northeast than 



that of the southwest. The detailed shell struct ures are also seen from the left panel, 



example, the "V-shape" knot at the southwest fiAschenbach fc Leahy 



1999 



Leahy 



200J). 



or 



From the left panel of flgure HI we found the flux distribution inside the Loop is far 
from what we expect in the uniform shell structure. This suggests the ambient density 
and the shell thickness varies considerably from region to region. Thus, we can study 
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the line-of-sight shell structure of the Loop. Considering the relation between the surface 
brightness and the plasma density, the flux of the low-kTf. component reflects the local 
density of the ISM. For example, the bright region in the northeast part is considered that 
the blast waves are expanding into the dense ISM there. In contrast, there is a low-flux 
region at the south of the Loop (see figure Hj). It suggests the ambi ent density there is 
extremely low compared with other areas of the Loop. As shown by lUchida et al.l (120081 ). 



we noticed that there is a large break in the south where the ISM density is very thin. In 
general, the velocity of the blast wave toward such tenuous ISM should become higher than 
other region. Therefore, it forms a blowout where the shell thickness must be thin. 

From figure HJ we also found a large low-flux region at slightly west of the Cygnus 
Loop center. Although our FOV does not cover the whole region, the structure is close to a 
circular form, and we estimated the diameter to be ~ 1.3°. The size is comparable to that of 
the south blowout. The existence of such large low-flux region suggests that it has a blowout 
str ucture along th e line of sig ht like the south blow out. This result confirms the prediction 



by 



Tsunemi et al. 



(120071 ) and 



Kimura et al. 



(120091 ) . From figure HI the northeast of the 
center also has lower flux than that of the surrounding region. It strongly indicates that 
the line-of-sight ambient density there is locally low as well as that in the south blowout. 
This region has a C-shape structure which could be explained by the superposition of the 
circular low-flux region and the bright region where the blast wave interacts with a small 
cloud. We estimate the diameter of this low-flux region to be ~ 30'. These results show the 



ambient density of the Cygnus Loop is quite different from region to region. 



4.3. Evidence of Cavity Explosion 

To put our result into perspective, we plotted the flux of each component as a function 
of radius R as shown in flgure [51 From this flgure, we found the flux of the high-fcTg 
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component (shown as crosses) decreases from the center to the outside, which reflects the 
spherical structure of the ejecta filled inside the Cygnus Loop. On the other hand, the 
flux of the low-A;Te component (circles) has a limb-brightening structure, as mentioned in 
the previous section. Furthermore the low-ZcTe flux at the southwest {R > 0) is totally 
lower than that at the northeast. While the high-fcTe flux distribution is approximately 
symmetric, the low-fcTe flux is a few times higher at the northeast than that at the 
southwest. In addition, looking at the inner region of the Loop, the flux distribution of the 
low-kTe component is declining from R = —50 to R = 50. This fact suggests the ambient 
density of the Cygnus Loop globally decreases from the northeast to the southwest. 

In order to estimate the ambient density more quantitatively, we calculated the EM 
of the low-fcTe component and plotted it as a function of R. Figure [6] shows the EM 
distribution of the low-fcTe component. We plotted the EM profiles from six rectangular 
regions with different azimuthal angles as shown in figure [6] (NE-A to NE-E and SW). That 
are shown in figure [3 We simulated the EM profile of the shell component derived from 
the Sedov solution with different ambient density no and estimated no by comparing our 
observations with the EM models. In this model, we assume the shock radius of the Cygnus 
Loop to be 13 pc and the ejecta is filled in 90% of it. The results are shown in figure [7] with 
red lines. We also show the best-fit models using the data only in the limb-brightening 
regions with green lines. As for the northeast regions, the EM profiles inside the Loop are 
close to the models of no=0.3-0.4 cm~^ (red) while the EM values at the limb-brightening 
regions are higher than these models. On the other hand, applying the data only in the 
limb-brightening regions (green), no increases to 0.7-0.9 cm~^. In any case, there are no 
Sedov models which agree with the EM profiles of the northeast part of the Cygnus Loop. 
The result is the same as the case of the southwest region while the ambient density uq is 
less than half of the northeast results. These results clearly show the Cygnus Loop can 



not be explained by a simple Sedov evolution model. 



McCray fc Snow 



(I1979I ) proposed 
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that the Cygnus Loop's S N explosion had occurred in a preexisting cavity and some ot 



studies also supported it (IHester fc Cox 



her 



1986 



Hester et al. 



1994 



Levenson et al. 



19971 ) 



Considering their results, it is natural that the EM distribution disagrees with a simple 
Sedov model, and thus, we concluded that our result also supports the cavity explosion 
as the origin of the Cygnus Loop from the standpoint of the X-ray spectral analysis. It 
should be noted that the Cygnus Loop is almost perfect circular in shape, although the EM 
(or flux) is globally higher in the northeast than that in the southwest. This fact strongly 
suggests that the northeast and the southwest blast waves should have hit the cavity wall 
very recently, and that the cavity-wall density is higher in the northeast than that in the 
southwest. 



5. Conclusion 

By analyzing the X-ray spectra, we clearly distinguished the ISM component from the 
ejecta component, and established a method to investigate the hne-of-sight shell structure. 
From the flux distribution of the ISM component, we found three low-fiux regions in the 
FOV; one is a well-known southwest blowout which is evidence of the cavity-wall break, 
and we also found other low-fiux regions at the west and the northeast of the Cygnus Loop 
center. From the EM distribution of the ISM component, we support that the Cygnus Loop 
is originated from a cavity explosion. Then, the ISM component, or cavity wall does not 
have an uniform structure but has a lot of breaks or tenuous regions. We also found that 
the condition of the surrounding cavity wall is not uniform; the density of it is globally 
higher in the northeast than that in the southwest. 



H.U. would like to thank Professor Jacco Vink and his students for many useful 
discussions and their hospitahty at Utrecht University. This work is partly supported by a 



- 15 - 



Grant-in- Aid for Scientific Research by the Ministry of Education, Culture, Sports, Science 
and Technology (16002004). H.U. and S.K. are supported by JSPS Research Fellowship for 
Young Scientists. 



-16- 





Fig. 1. — Top: ROSAT HRI image of the entire Cygnus Loop. The circles and rectangles 
represent our FOV of the XMM-Newton MOS and the Suzaku XIS, respectively. Bottom: 
Same as the left panel, but for overlaid with the spectral extraction regions with small 
rectangles. 
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Fig. 2. — Example XISl spectra from the regions where the flux of the swept-up matter 
is high (region-A: left two panels) and low (region-B: right two panels), respectively (see 
figur^I]). The best-fit curves for the single-component VNEI models are shown by solid 
black lines in the top two panels. Bottom two panels are the same as the top panels, but for 
the fitting results with the two-component VNEI models. In the bottom panels, blue and 
red lines represent the high-ZcTg component and the low-ZcTg component, respectively. The 
residuals are shown in lower panels. 
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Fig. 3. — Our FOV and the electron temperature distribution of the \ow-kTe component 
overlaid with the white contour from the ROSAT HRI image. The images are smoothed by 
Gaussian kernel of a = 2.8'. The values are in units of keV. 
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0.002 0.006 0.01 



Fig. 4. — 0.2-3.0 keV flux distribution of tlie low-kT^ (left) and tlie higfi-fcTe (right) compo- 
nent in logaritlimic scales overlaid with the white contour of the ROSAT HRI image. The 
images are smoothed by Gaussian kernel of a = 2.8'. The values are in units of counts 
cm~^s~^arcmin~^ and the scale parameters correspond with each other. Blue and red corre- 
spond to ~ 10"'' and ~ 10~^ counts cm^^s^^arcmin"^, respectively. 
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Radius [arcmin] 



Fig. 5. — Averaged flux profile as a function of R. The circles and triangles represent the 
flux of low-A;Te and high-fcTg components, respectively. 



- 21 - 




5E+20 1E+21 2E+21 



Fig. 6. — EM distribution of the low-kT^ component in logarithmic scales overlaid with the 
white contour of the ROSAT HRI image. 
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Fig. 7. — EM profiles as a function of R calculated from the data in the rectangular regions 
shown in figure [6l The EM profiles based on the Sedov model and the estimated ambient 
densities no are shown in red and green (see text). 
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Table 1. Summary of the 41 observations 



Obs. ID Obs. Date RA, DEC (J2000) Position Angle Effective Exposure 



Suzaku Observations 



501012010 


(PI) 


2007-11- 


■14 


20''54"^07.6^ 


31°57'22.0" 


240°.0 


9.8 ksec 


501013010 


(P2) 


2007-11- 


-14 


20^53°'08.5^ 


3r45'40.3" 


240°. 


16.4 kscc 


501014010 


(P3) 


2007-11- 


-14 


20^52°^09.9^ 


31°36'43.4" 


240°.0 


16.9 ksec 


501015010 


(P4) 


2007-11- 


■14 


20''51"'11.8^ 


31°22'08.4" 


240°.0 


18.3 ksec 


501016010 


(P5) 


2007-11- 


-15 


20''50'"11.3^ 


31°10'48.0" 


240°. 


19.3 ksec 


501017010 


(P6) 


2007-11- 


-11 


20M9'"11.3^ 


30°59'27.6" 


240°. 


28.7 ksec 


501018010 


(P7) 


2007-11- 


■12 


20M8"'18.7^ 


30°46'33.6" 


240°.0 


21.0 ksec 


501028010 


(P8) 


2006-05- 


■13 


20''55"'56.3^ 


31°28'56.2" 


62°.5 


4.9 ksec 


501019010 


(P9) 


2007-11- 


-12 


20M7'^14.2% 


30°36'10.8" 


240°. 


16.2 ksec 


501020010 


(PIO) 


2007-11- 


-13 


20M6"'20.8^ 


30°23'22.6" 


240°. 


14.7 ksec 


503055010 


(Pll) 


2008-05- 


■09 


20M9"'48.7^ 


31°30'18.0" 


50°.0 


22.2 ksec 


501029010 


(P12) 


2006-05- 


■09 


20''55"'00.0^ 


31°15'46.8" 


62°. 1 


13.2 ksec 


501030010 


(P13) 


2006-05- 


-10 


20^53°^59.3^ 


31°03'39.6" 


68°. 2 


13.9 ksec 


501031010 


(P14) 


2006-05- 


■12 


20^52"^58.8^ 


30°51'32.4" 


62°.4 


18.2 ksec 


501032010 


(P15) 


2006-05- 


■25 


20''51"'58.6^ 


30°39'10.8" 


62°.0 


17.4 ksec 


501033010 


(P16) 


2006-05- 


-22 


20^^50°'58.8% 


30°27'00.0" 


62°. 


20.0 ksec 


501034010 


(P17) 


2006-05- 


-22 


20M8"^49.7^ 


30°00'21.6" 


62°.0 


13.9 ksec 


501035010 


(P18) 


2006-12- 


■18 


20M8'"16.2^ 


29°42'07.2" 


237°.5 


11.2 ksec 


501036010 


(P19) 


2006-12- 


■18 


2Qh47m^73s^ 


30°04'21.4" 


237°. 5 


11.8 ksec 


503056010 


(P20) 


2008-05- 


-10 


20M8"'00.0^ 


31°10'30.0" 


50°. 


22.5 ksec 



Table 1 — Continued 



Obs. ID Obs. Date RA, DEC (J2000) Position Angle Effective Exposure 



503057010 (P21) 


2008-06- 


02 


20^ 


52"^43.8^ 


32°26'19.0" 


61°.9 


16.2 


ksec 


503058010 (P22) 


2008-06- 


03 


20^ 


51^^17.2% 


32°25'24.6" 


61°.4 


19.3 


ksec 


503059010 (P23) 


2008-06- 


03 


20*^ 


49"'50.6^ 


32°21'50.8" 


61°.9 


19.5 


ksec 


503060010 (P24) 


2008-06- 


04 


20*^ 


48"'28.2^ 


32°ir44.5" 


61°.4 


18.5 


ksec 


503061010 (P25) 


2008-06- 


04 


20^ 


47^22.7^, 


32°10'22.8" 


60°.9 


26.0 


ksec 


503062010 (P26) 


2008-05- 


13 


20^ 


56'^26.5^ 


30°19'55.2" 


49°.8 


16.9 


ksec 


503063010 (P27) 


2008-05- 


13 


20*^ 


55'"16.3^ 


30°01'44.0" 


49°.6 


22.8 


ksec 


503064010 (P28) 


2008-05- 


14 


20^ 


53"'51.6^ 


29°54'42.5" 


49°. 1 


18.2 


ksec 


500020010 (NEl) 


2005-11- 


23 


20^" 


56'^48.9% 


3r56'54.8" 


223°.0 


20.4 


ksec 


500021010 (NE2) 


2005-11- 


24 


20*^ 


55"'56.0^ 


31°56'53.2" 


223°.0 


21.4 


ksec 


500022010 (NE3) 


2005-11- 


29 


20*^ 


55"^05.6^ 


32°10'35.4" 


222°.9 


21.7 


ksec 


500023010 (NE4) 


2005-11- 


30 


20^ 


54"^03.8^ 


32°21'47.9" 


221°.2 


25.3 


ksec 



XMM-Newton Observations 



0082540101 (Pos-1) 


2002-11- 


25 


20^^55^^23.6% 


31°46'17.0" 


241°. 7 


14.7 ksec 


0082540201 (Pos-2) 


2002-12- 


03 


20'^54"^07.2% 


31°30'51.1" 


241°.7 


14.4 ksec 


0082540301 (Pos-3) 


2002-12- 


05 


20''52"^51.1% 


31°15'25.7" 


241°.7 


11.6 ksec 


0082540401 (Pos-4) 


2002-12- 


07 


20^51"'34.7% 


31°00'00.0" 


241°.7 


4.9 ksec 


0082540501 (Pos-5) 


2002-12- 


09 


20^50^^18.4% 


30°44'34.3" 


231°.4 


12.6 ksec 


0082540601 (Pos-6) 


2002-12- 


11 


20M9"'02.0% 


30°29'08.6" 


241°.7 


11.5 ksec 


0082540701 (Pos-7) 


2002-12- 


13 


20M7"'45.8% 


30°13'42.9" 


241°.7 


13.7 ksec 


0405490101 (Pos-8) 


2006-05- 


13 


20''50°^32.2% 


30°11'00.0" 


69°.9 


6.5 ksec 
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Table 1 — Continued 



Obs. ID Obs. Date RA, DEC (J2000) Position Angle Effective Exposure 



0405490201 (Pos-9) 2006-05-13 20M9"^54.2^ 29°42'25.0" 69°.8 



3.6 ksec 
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Table 2: Spectral fit parameters 



single- component VNEI model two-component VNEI model 





region A 


region B 


region A 


region B 


Nh [lO^Ocm-^] 


1.8 ± 0.3 


3.4 ± 0.3 


5.2 ± 0.2 


7.0 ± 0.3 








Low-kTe component: 


kT^ [keV] 


0.59 ± 0.03 


0.42 ± 0.02 


0.24 ± 0.01 


0.12 ± 0.01 


C 


0.24 ± 0.05 


0.96 ± 0.21 


0.27 


(flxed) 


N 


0.22 ± 0.05 


0.09 ± 0.03 


0.10 


(fixed) 





0.23 ± 0.02 


0.14 ± 0.02 


0.11 


(fixed) 


Ne 


0.44 ± 0.04 


0.31 ± 0.03 


0.21 


(flxed) 


Mg 


0.26 ± 0.03 


0.24 ± 0.03 


0.17 


(flxed) 


Si 


0.27 ± 0.06 


0.30 ± 0.06 


0.34 


(flxed) 


S 


(=Si) 


(=Si) 


0.17 


(flxed) 


Fe(=Ni) 


0.36 ± 0.04 


0.22 ± 0.02 


0.20 


(flxed) 


log r 


10.42 ± 0.03 


10.82 1^8^ 


11.32 t'oil 


< 12 


flux [counts cm~^s~^arcmin~^] 


8.90 X 10"^ 


4.34 X 10"^ 


7.49 X 10"^ 
High-kTe 


2.18 X 10"^ 
component: 


kTe [keV] 






0.88 ± 0.13 


0.43 ± 0.02 


0(=C=N) 






0.34 ± 0.13 


0.38 ± 0.07 


Ne 






0.82 ± 0.26 


0.74 ± 0.12 


Mg 






0.56 ± 0.19 


0.55 ± 0.10 


Si(=S) 






1.28 ± 0.42 


0.79 ± 0.16 


Fe(=Ni) 






< 1.43 


0.48 ± 0.08 


log T 






10.66 


11.11 ± 0.06 


flux [counts cm~^s~^arcmin~^] 






1.41 X 10-^ 


2.16 X 10-^ 


xVdof 


1043/739 


728/548 


868/738 


637/547 
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